Longitudinal genetic dynamics of weaning index and implications for cow-calf production efficiency.

In beef cattle, the selection for higher weights at young ages has been questioned with the argument that this criterion may increase the adult weight of cows, resulting in higher costs. Therefore, selection criteria should be employed to increase weights at young ages with minimal impact on the adult weight of cows. Additionally, the relationship between measures of cow production efficiency and other well-established selection criteria in breeding programs remains poorly understood. The objective of this study was to longitudinally evaluate the relationship between the weaning index (WIndex) as a measure of efficiency and growth traits of the cows. Possible changes over time in WIndex due to selection applied for yearling weight (YW) were also investigated. The WIndex was proposed to maximize genetic response in the weaning weight of the calf while maintaining genetic gain in BW of the cow at zero. A random regression model was adopted to estimate correlations between WIndex, BW, hip height (HH), and body condition score (BCS) using records of Nelore cows from three lines. Genetic trends were calculated for the control line (NeC) and lines selected for greater YW (NeS and NeT). The age of 3 years was the most critical for the weaning efficiency of the cows. At this stage, young cows are still growing and wean lighter calves than their adult counterparts. The genetic correlation estimates between WIndex and BW (-0.58 to 0.04), HH (-0.05 to -0.34), and BCS (-0.51 to -0.17) were close to zero or negative. BW and HH were strongly correlated genetically across all ages (0.73-0.76). In general, HH exhibited a weak and negative genetic relationship with BCS. The genetic correlation between BW and BCS was stronger for advanced ages (0.45-0.68). In lines selected for YW, important increases in WIndex were observed. However, NeS has been selected since the 1980s until the present for YW, and thus, it showed a more pronounced trend of increasing BW and, consequently, a more modest trend of increasing WIndex compared to NeT. In contrast, WIndex exhibited a trend close to zero for NeC. In this context, monitoring HH and BCS can be useful to avoid losses in the weaning efficiency of cows. Furthermore, we suggest that one way to mitigate efficiency losses in calf production could involve stabilizing the BW of cows and increasing the weaning weight of calves using the WIndex.

Genetic component of sensitivity to heat stress for nonreturn rate of Brazilian Holstein cattle.

The objectives of the present study were: 1) to investigate variation in the genetic component of heat stress for nonreturn rate at 56 days after first artificial insemination (NR56); 2) to identify and characterize the genotype by environment interaction (G × E) due to heat stress for NR56 of Brazilian Holstein cattle. A linear random regression model (reaction norm model) was applied to 51,748 NR56 records of 28,595 heifers and multiparous cows. The decline in NR56 due to heat stress was more pronounced in milking cows compared to heifers. The age of females at first artificial insemination and temperature-humidity index (THI) exerted an important influence on the genetic parameters of NR56. Several evidence of G × E on NR56 were found as the high slope/intercept ratio and frequent intersection of reaction norms. Additionally, the genetic correlation between NR56 at opposite extremes of the THI scale reached estimates below zero, indicating that few of the same genes are responsible for NR56 under conditions of thermoneutrality and heat stress. The genetic evaluation and selection for NR56 in Holstein cattle reared under (sub)tropical conditions should therefore take into consideration the genetic variation on age at insemination and G × E due to heat stress.

Genetics of heat tolerance for milk yield and quality in Holsteins.

Tropical and sub-tropical climates are characterized by high temperature and humidity, during at least part of the year. Consequently, heat stress is common in Holstein cattle and productive and reproductive losses are frequent. Our objectives were as follows: (1) to quantify losses in production and quality of milk due to heat stress; (2) to estimate genetic correlations within and between milk yield (MY) and milk quality traits; and (3) to evaluate the trends of genetic components of tolerance to heat stress in multiple lactations of Brazilian Holstein cows. Thus, nine analyses using two-trait random regression animal models were carried out to estimate variance components and genetic parameters over temperature-humidity index (THI) values for MY and milk quality traits (three lactations: MY×fat percentage (F%), MY×protein percentage (P%) and MY×somatic cell score (SCS)) of Brazilian Holstein cattle. It was demonstrated that the effects of heat stress can be harmful for traits related to milk production and milk quality of Holstein cattle even though most herds were maintained in a modified environment, for example, with fans and sprinklers. For MY, the effect of heat stress was more detrimental in advanced lactations (-0.22 to -0.52 kg/day per increase of 1 THI unit). In general, the mean heritability estimates were higher for lower THI values and longer days in milk for all traits. In contrast, the heritability estimates for SCS increased with increasing THI values in the second and third lactation. For each trait studied, lower genetic correlations (different from unity) were observed between opposite extremes of THI (THI 47 v. THI 80) and in advanced lactations. The genetic correlations between MY and milk quality trait varied across the THI scale and lactations. The genotype×environment interaction due to heat stress was more important for MY and SCS, particularly in advanced lactations, and can affect the genetic relationship between MY and milk quality traits. Selection for higher MY, F% or P% may result in a poor response of the animals to heat stress, as a genetic antagonism was observed between the general production level and specific ability to respond to heat stress for these traits. Genetic trends confirm the adverse responses in the genetic components of heat stress over the years for milk production and quality. Consequently, the selection of Holstein cattle raised in modified environments in both tropical and sub-tropical regions should take into consideration the genetic variation in heat stress.

Inbreeding depression in Zebu cattle traits.

The productivity of herds may be negatively affected by inbreeding depression, and it is important to know how intense is this effect on the livestock performance. We performed a comprehensive analysis involving five Zebu breeds reared in Brazil to estimate inbreeding depression in productive and reproductive traits. Inbreeding depression was estimated for 13 traits by including the individual inbreeding rate as a linear covariate in the standard genetic evaluation models. For all breeds and for almost all traits (no effect was observed on gestation length), the performance of the animals was compromised by an increase in inbreeding. The average inbreeding depression was -0.222% and -0.859% per 1% of inbreeding for linear regression coefficients scaled on the percentage of mean (ßm ) and standard deviation (ßσ ), respectively. The means for ßm (and ßσ ) were -0.269% (-1.202%) for weight/growth traits and -0.174% (-0.546%) for reproductive traits. Hence, inbreeding depression is more pronounced in weight/growth traits than in reproductive traits. These findings highlight the need for the management of inbreeding in the respective breeding programmes of the breeds studied here.

Genetic variation of the weaning weight of beef cattle as a function of accumulated heat stress.

The objective of this study was to identify the genetic variation in the weaning weight (WW) of beef cattle as a function of heat stress. The WWs were recorded at approximately 205 days of age in three Brazilian beef cattle populations: Nelore (93,616), Brangus (18,906) and Tropical Composite (62,679). In view of the cumulative nature of WW, the effect of heat stress was considered as the accumulation of temperature and humidity index units (ACTHI) from the animal's birth to weaning. A reaction norm model was used to estimate the (co)variance components of WW across the ACTHI scale. The accumulation of THI units from birth to weaning negatively affected the WW. The definition of accumulated THI units as an environmental descriptor permitted to identify important genetic variation in the WW as a function of heat stress. As evidence of genotype by environment interaction, substantial heterogeneity was observed in the (co)variance components for WW across the environmental gradient. In this respect, the best animals in less stressful environments are not necessarily the best animals in more stressful environments. Furthermore, the response to selection for WW is expected to be lower in more stressful environments.

Random regression analysis of test-day milk yields in the first and second lactations of Brazilian Gyr cows.

The objective of the present study was to estimate the genetic parameters for test-day milk yields (TDMY) in the first and second lactations using random regression models (RRM) in order to contribute to the application of these models in genetic evaluation of milk yield in Gyr cattle. A total of 53,328 TDMY records from 7118 lactations of 5853 Gyr cows were analyzed. The model included the direct additive, permanent environmental, and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cows at calving were included as fixed effects. A random regression model fitting fourth-order Legendre polynomials for additive genetic and permanent environmental effects, with five classes of residual variance, was applied. In the first lactation, the heritabilities increased from early lactation (0.26) until TDMY3 (0.38), followed by a decrease until the end of lactation. In the second lactation, the estimates increased from the first (0.29) to the fifth test day (0.36), with a slight decrease thereafter, and again increased on the last two test days (0.34 and 0.41). There were positive and high genetic correlations estimated between first-lactation TDMY and the remaining TDMY of the two lactations. The moderate heritability estimates, as well as the high genetic correlations between half the first-lactation TDMY and all TDMY of the two lactations, suggest that the selection based only on first lactation TDMY is the best selection strategy to increase milk production across first and second lactations of Gyr cows.

Genetic parameters of test-day milk yield in Guzerá cattle under tropical conditions.

The objective of this study was to estimate genetic parameters for test-day milk yield (TDMY) in Guzerá cows using random regression models. Additive and permanent environmental random effects were modeled by random regression on fourth- and fifth-order orthogonal Legendre polynomials, respectively. The residual variances were heterogeneous, with seven classes. Heritability estimates for TDMY ranged from 0.24 to 0.52, with higher heritabilities for yields during early lactation. Genetic correlations between TDMYs ranged from -0.03 to 0.95. The phenotypic and permanent environmental correlations were all positive, and the highest estimates were between adjacent TDMYs. The results suggest that TDMYs obtained with random regression models may be used as selection criteria for Guzerá cattle.

Multivariate analysis of test-day and total milk yield in goats.

The objective of this study was to estimate genetic parameters for 305-day cumulative milk yield (MY305) and its association with test-day milk yield (TDMY) in Saanen and Alpine goats in order to provide information that allows the use of TDMY as selection criteria. This was done using standard multi-trait and reduced rank models. Data from 1157 lactations, including the first three kiddings, and 5435 test-day records from 683 Saanen and 449 Alpine goats were used. MY305 was analyzed together with TDMY by multi-trait analysis, from the first to tenth test-day, using records of the first three lactations as repeated measures. Three multi-trait models were used: a standard (SM) and two reduced rank models that fitted the first two (PC2) and three (PC3) genetic principal components. Akaike and Schwarz Bayesian information criteria were used to compare models. Heritability for TDMY estimated with the SM ranged from 0.20 to 0.66, whereas the range calculated from the PC2 model was 0.16 to 0.63. Genetic correlations between TDMY and MY305 were positive and moderate to high, ranging from 0.56 to 0.98 when estimated with the SM, and 0.91 to 1.00 when estimated with the PC2. The standard multi-trait model produced estimates that were more accurate than the reduced rank models. Although the SM provided the worst fit according to the two model selection criteria, it was the best in this dataset.

Detrimental effect of selection for milk yield on genetic tolerance to heat stress in purebred Zebu cattle: Genetic parameters and trends.

In an attempt to determine the possible detrimental effects of continuous selection for milk yield on the genetic tolerance of Zebu cattle to heat stress, genetic parameters and trends of the response to heat stress for 86,950 test-day (TD) milk yield records from 14,670 first lactations of purebred dairy Gir cows were estimated. A random regression model with regression on days in milk (DIM) and temperature-humidity index (THI) values was applied to the data. The most detrimental effect of THI on milk yield was observed in the stage of lactation with higher milk production, DIM 61 to 120 (-0.099kg/d per THI). Although modest variations were observed for the THI scale, a reduction in additive genetic variance as well as in permanent environmental and residual variance was observed with increasing THI values. The heritability estimates showed a slight increase with increasing THI values for any DIM. The correlations between additive genetic effects across the THI scale showed that, for most of the THI values, genotype by environment interactions due to heat stress were less important for the ranking of bulls. However, for extreme THI values, this type of genotype by environment interaction may lead to an important error in selection. As a result of the selection for milk yield practiced in the dairy Gir population for 3 decades, the genetic trend of cumulative milk yield was significantly positive for production in both high (51.81kg/yr) and low THI values (78.48kg/yr). However, the difference between the breeding values of animals at high and low THI may be considered alarming (355kg in 2011). The genetic trends observed for the regression coefficients related to general production level (intercept of the reaction norm) and specific ability to respond to heat stress (slope of the reaction norm) indicate that the dairy Gir population is heading toward a higher production level at the expense of lower tolerance to heat stress. These trends reflect the genetic antagonism between production and tolerance to heat stress demonstrated by the negative genetic correlation between these components (-0.23). Monitoring trends of the genetic component of heat stress would be a reasonable measure to avoid deterioration in one of the main traits of Zebu cattle (i.e., high tolerance to heat stress). On the basis of current genetic trends, the need for future genetic evaluation of dairy Zebu animals for tolerance to heat stress cannot be ruled out.

Reduced-rank models of growth and reproductive traits in Nelore cattle.

In beef cattle genetic evaluation, principal component models of the additive genetic effect could be used to incorporate several traits in the same analysis, without an important increase in the number of parameters to be estimated. In this study, multitrait (MT) and reduced-rank models were compared for their ability to estimate parameters and predict breeding values for weaning weight, yearling weight, weaning hip height, yearling hip height, weaning to yearling weight gain, scrotal circumference, and age at the first calving. Data obtained were from 74,388 Nelore animals, born to 1441 sires and 28,502 cows. Six analyses were performed using a MT model that incorporated all the traits simultaneously and five reduced-rank models for the genetic additive direct (co)variance matrix, fitting the first one (PC1), two (PC2), three (PC3), four (PC4), and five (PC5) principal components. The model considering the first three principal components (PC3) provided the best fit. Direct and maternal heritability and the respective standard errors obtained from the MT and PC3 models were similar. In general, the PC3 model provided slightly stronger genetic correlations between the traits when compared with those obtained with the MT model. The rank correlations between the breeding values predicted with the MT and PC3 models for the different traits ranged from 0.93 to 0.99. When 2% and 10% of the best sires were selected on the basis of breeding values predicted by the MT model, the degree of concordance with the PC3 model ranged from 86% to 97%. The first three principal components explained most of the genetic variation among animals, suggesting that major changes should not be expected in the sire's classification on the basis of breeding values predicted for growth and reproductive traits. Models of principal components could be used for beef cattle genetic evaluation, especially when considering several economic traits in the same analysis.

Two-trait random regression model to estimate the genetic association of scrotal circumference with female reproductive performance in Nelore cattle.

In an attempt to determine when scrotal circumference (SC) could be a reasonable indicator of female reproductive performance, a series of two-trait random regression model (regression for SC on age at measurement) using Gibbs sampling was applied to field data of Nelore cattle raised in a tropical environment. The female traits evaluated were age at first calving (AFC), first calving interval (FCI), heifer pregnancy (HP), and stayability (STAY). The posterior means of heritability of female traits ranged from 0.15 for AFC to 0.46 for HP and were about 0.50 for SC. The posterior means of genetic correlations between SC and AFC, FCI, HP, STAY were up to -0.70, -0.25, 0.48, and 0.29, respectively. Genetically, SC could be a reasonable indicator of female puberty (e.g., HP) as long as it is measured at a young age (400-440 days). However, for female reproductive traits measured at an older age (e.g., STAY), SC is merely a modest or even poor indicator. The use of sire expected progeny differences for female reproductive traits will be more effective than the use of expected progeny differences for SC to improve the reproductive performance of female cattle.

Random regression models using different functions to model milk flow in dairy cows.

We analyzed 75,555 test-day milk flow records from 2175 primiparous Holstein cows that calved between 1997 and 2005. Milk flow was obtained by dividing the mean milk yield (kg) of the 3 daily milking by the total milking time (min) and was expressed as kg/min. Milk flow was grouped into 43 weekly classes. The analyses were performed using a single-trait Random Regression Models that included direct additive genetic, permanent environmental, and residual random effects. In addition, the contemporary group and linear and quadratic effects of cow age at calving were included as fixed effects. Fourth-order orthogonal Legendre polynomial of days in milk was used to model the mean trend in milk flow. The additive genetic and permanent environmental covariance functions were estimated using random regression Legendre polynomials and B-spline functions of days in milk. The model using a third-order Legendre polynomial for additive genetic effects and a sixth-order polynomial for permanent environmental effects, which contained 7 residual classes, proved to be the most adequate to describe variations in milk flow, and was also the most parsimonious. The heritability in milk flow estimated by the most parsimonious model was of moderate to high magnitude.

Genetic relationships among traits related to reproduction and growth of Nelore females.

The objective of the study presented here was to analyze the genetic relationships among heifer pregnancy (HP), age at first calving (AFC), stayability (STAY), average annual productivity of the cow, in kilograms of weaned calf per cow per year (PRODAM), postweaning weight gain (PWG), and hip height (HH) of Nelore females from 12 Brazilian herds. (Co)variance components were obtained by six-trait animal model using Gibbs sampling. The posterior mean of the heritability estimates were 0.37, 0.18, 0.19, 0.16, 0.21, and 0.37 for HP, AFC, STAY, PRODAM, PWG, and HH, respectively. In general, the genetic correlations were strong between traits related to reproduction, for example, -0.85 between HP and AFC, and 0.94 between STAY and PRODAM. Weak genetic correlations were obtained between reproductive and growth traits (absolute values ranging from 0.02 to 0.30). Although weak, the genetic correlations between PWG and reproductive traits were favorable, whereas the genetic correlations between HH and reproductive traits were close to zero and slightly unfavorable for HP, AFC, and STAY. An increase of HH is therefore expected to have little or no negative effect on the reproductive performance of females. The posterior mean of genetic correlation between PWG and HH was moderate (0.50). On the basis of the heritability, genetic correlation estimates, and time to obtain data, HP and PRODAM seems to show the best potential as selection criteria to improve the productive and reproductive performance of Nelore females. In principle, it is possible to select for increased PWG without compromising the reproduction of Nelore females. However, selection for PWG may result in an increase of female HH as a correlated response, a fact that could increase management costs in advanced generations of selection. In the light of the results, all traits studied here can be used as selection criteria and there is no strong evidence of genetic antagonism among traits related to reproduction and growth of Nelore females.

Estimates of genetic parameters for total milk yield over multiple ages in Brazilian Murrah buffaloes using different models.

The objective of this study was to estimate variance components and genetic parameters for accumulated 305-day milk yield (MY305) over multiple ages, from 24 to 120 months of age, applying random regression (RRM), repeatability (REP) and multi-trait (MT) models. A total of 4472 lactation records from 1882 buffaloes of the Murrah breed were utilized. The contemporary group (herd-year-calving season) and number of milkings (two levels) were considered as fixed effects in all models. For REP and RRM, additive genetic, permanent environmental and residual effects were included as random effects. MT considered the same random effects as did REP and RRM with the exception of permanent environmental effect. Residual variances were modeled by a step function with 1, 4, and 6 classes. The heritabilities estimated with RRM increased with age, ranging from 0.19 to 0.34, and were slightly higher than that obtained with the REP model. For the MT model, heritability estimates ranged from 0.20 (37 months of age) to 0.32 (94 months of age). The genetic correlation estimates for MY305 obtained by RRM (L23.res4) and MT models were very similar, and varied from 0.77 to 0.99 and from 0.77 to 0.99, respectively. The rank correlation between breeding values for MY305 at different ages predicted by REP, MT, and RRM were high. It seems that a linear and quadratic Legendre polynomial to model the additive genetic and animal permanent environmental effects, respectively, may be sufficient to explain more parsimoniously the changes in MY305 genetic variation with age.

Genotype by production environment interaction for birth and weaning weights in a population of composite beef cattle.

The objectives of the present study were: (1) to evaluate the importance of genotype × production environment interaction for the genetic evaluation of birth weight (BW) and weaning weight (WW) in a population of composite beef cattle in Brazil, and (2) to investigate the importance of sire × contemporary group interaction (S × CG) to model G × E and improve the accuracy of prediction in routine genetic evaluations of this population. Analyses were performed with one, two (favorable and unfavorable) or three (favorable, intermediate, unfavorable) different definitions of production environments. Thus, BW and WW records of animals in a favorable environment were assigned to either trait 1, in an intermediate environment to trait 2 or in an unfavorable environment to trait 3. The (co)variance components were estimated using Gibbs sampling in single-, bi- or three-trait animal models according to the definition of number of production environments. In general, the estimates of genetic parameters for BW and WW were similar between environments. The additive genetic correlations between production environments were close to unity for BW; however, when examining the highest posterior density intervals, the correlation between favorable and unfavorable environments reached a value of only 0.70, a fact that may lead to changes in the ranking of sires across environments. The posterior mean genetic correlation between direct effects was 0.63 in favorable and unfavorable environments for WW. When S × CG was included in two- or three-trait analyses, all direct genetic correlations were close to unity, suggesting that there was no evidence of a genotype × production environment interaction. Furthermore, the model including S × CG contributed to prevent overestimation of the accuracy of breeding values of sires, provided a lower error of prediction for both direct and maternal breeding values, lower squared bias, residual variance and deviance information criterion than the model omitting S × CG. Thus, the model that included S × CG can therefore be considered the best model on the basis of these criteria. The genotype × production environment interaction should not be neglected in the genetic evaluation of BW and WW in the present population of beef cattle. The inclusion of S × CG in the model is a feasible and plausible alternative to model the effects of G × E in the genetic evaluations.

Genotype by environment interaction and model comparison for growth traits of Santa Ines sheep.

The objectives of the present study were to compare alternative models for the genetic evaluation and assess the importance of genotype by environment interaction (G×E) in the estimation of genetic parameters and genetic evaluation of birth weight (BW), weight at 60 days of age (W60) and weight at 180 days of age (W180) of Santa Ines sheep. Data comprise 7622 BW, 4673 W60 and 2830 W180 records from animals born in 44 Brazilian herds. Four models were used for the analyses: animal model (AM) with homogeneous residual variance (1), or heterogeneous residual variance (2), hierarchical reaction norms model (HRNM) with homogeneous (1) or heterogeneous residual variance (2). The models that best fit the BW, W60 and W180 data were AM2, HRNM1 and HRNM2 respectively. Thus, models for genetic evaluation that consider heterogeneity of variances are recommended to evaluate growth traits of sheep. The correlation between intercept and slope of the HRNM was higher than 0.70 for all traits studied, indicating that animals with higher average breeding values responded better to improvement in environmental conditions, a fact characterizing the scale effect of G×E. Therefore, G×E is an important factor to be considered in the estimation of genetic parameters and genetic evaluation of growth traits of sheep.

Genetic associations among average annual productivity, growth traits, and stayability: a parallel between Nelore and composite beef cattle.

This study was conducted to examine the relationship among average annual productivity of the cow (PRODAM), yearling weight (YW), postweaning BW gain (PWG), scrotal circumference (SC), and stayability in the herd for at least 6 yr (STAY) of Nelore and composite beef cattle. Measurements were taken on animals born between 1980 and 2010 on 70 farms located in 7 Brazilian states. Estimates of heritability and genetic and environmental correlations were obtained by Bayesian approach with 5-trait animal models. Genetic trends were estimated by regressing means of estimated breeding values by year of birth. The heritability estimates were between 0.14 and 0.47. Estimates of genetic correlation among female traits (PRODAM and STAY) and growth traits ranged from -0.02 to 0.30. Estimates of genetic correlations ranged from 0.23 to 0.94 among growth traits indicating that selection for these traits could be successful in tropical breeding programs. Genetic correlations among all traits were favorable and simultaneous selection for growth, productivity, and stayability is therefore possible. Genetic correlation between PRODAM and STAY was 0.99 and 0.85 for Nelore and composite cattle, respectively. Therefore, PRODAM and STAY might be influenced by many of the same genes. The inclusion of PRODAM instead of STAY as a selection criterion seems to be more advantageous for tropical breeding programs because the generation interval required to obtain accurate estimates of genetic merit for PRODAM is shorter. Average annual genetic changes were greater in Nelore than in composite cattle. This was not unexpected because the breeding program of composite cattle included a large number of farms, different production environments, and genetic level of the herds and breeds. Thus, the selection process has become more difficult in this population.

Random regression models using Legendre polynomials or linear splines for test-day milk yield of dairy Gyr (Bos indicus) cattle.

Studies investigating the use of random regression models for genetic evaluation of milk production in Zebu cattle are scarce. In this study, 59,744 test-day milk yield records from 7,810 first lactations of purebred dairy Gyr (Bos indicus) and crossbred (dairy Gyr × Holstein) cows were used to compare random regression models in which additive genetic and permanent environmental effects were modeled using orthogonal Legendre polynomials or linear spline functions. Residual variances were modeled considering 1, 5, or 10 classes of days in milk. Five classes fitted the changes in residual variances over the lactation adequately and were used for model comparison. The model that fitted linear spline functions with 6 knots provided the lowest sum of residual variances across lactation. On the other hand, according to the deviance information criterion (DIC) and bayesian information criterion (BIC), a model using third-order and fourth-order Legendre polynomials for additive genetic and permanent environmental effects, respectively, provided the best fit. However, the high rank correlation (0.998) between this model and that applying third-order Legendre polynomials for additive genetic and permanent environmental effects, indicates that, in practice, the same bulls would be selected by both models. The last model, which is less parameterized, is a parsimonious option for fitting dairy Gyr breed test-day milk yield records.

Genetic parameters for first lactation test-day milk flow in Holstein cows.

Genetic parameters for test-day milk flow (TDMF) of 2175 first lactations of Holstein cows were estimated using multiple-trait and repeatability models. The models included the direct additive genetic effect as a random effect and contemporary group (defined as the year and month of test) and age of cow at calving (linear and quadratic effect) as fixed effects. For the repeatability model, in addition to the effects cited, the permanent environmental effect of the animal was also included as a random effect. Variance components were estimated using the restricted maximum likelihood method in single- and multiple-trait and repeatability analyses. The heritability estimates for TDMF ranged from 0.23 (TDMF 6) to 0.32 (TDMF 2 and TDMF 4) in single-trait analysis and from 0.28 (TDMF 7 and TDMF 10) to 0.37 (TDMF 4) in multiple-trait analysis. In general, higher heritabilities were observed at the beginning of lactation until the fourth month. Heritability estimated with the repeatability model was 0.27 and the coefficient of repeatability for first lactation TDMF was 0.66. The genetic correlations were positive and ranged from 0.72 (TDMF 1 and 10) to 0.97 (TDMF 4 and 5). The results indicate that milk flow should respond satisfactorily to selection, promoting rapid genetic gains because the estimated heritabilities were moderate to high. Higher genetic gains might be obtained if selection was performed in the TDMF 4. Both the repeatability model and the multiple-trait model are adequate for the genetic evaluation of animals in terms of milk flow, but the latter provides more accurate estimates of breeding values.

Short communication: Principal components and factor analytic models for test-day milk yield in Brazilian Holstein cattle.

A total of 46,089 individual monthly test-day (TD) milk yields (10 test-days), from 7,331 complete first lactations of Holstein cattle were analyzed. A standard multivariate analysis (MV), reduced rank analyses fitting the first 2, 3, and 4 genetic principal components (PC2, PC3, PC4), and analyses that fitted a factor analytic structure considering 2, 3, and 4 factors (FAS2, FAS3, FAS4), were carried out. The models included the random animal genetic effect and fixed effects of the contemporary groups (herd-year-month of test-day), age of cow (linear and quadratic effects), and days in milk (linear effect). The residual covariance matrix was assumed to have full rank. Moreover, 2 random regression models were applied. Variance components were estimated by restricted maximum likelihood method. The heritability estimates ranged from 0.11 to 0.24. The genetic correlation estimates between TD obtained with the PC2 model were higher than those obtained with the MV model, especially on adjacent test-days at the end of lactation close to unity. The results indicate that for the data considered in this study, only 2 principal components are required to summarize the bulk of genetic variation among the 10 traits.